We present the first reliable calculation of the collec- 
tive mode structure of a strongly coupled electronic bilayer. 
The calculation is based on a classical model through the 
3 rd frequency-moment-sum-rule preserving Quasi Localized 
Charge Approximation, using the recently calculated Hyper- 
netted Chain pair correlation functions. The spectrum shows 
an energy gap at k = and the absence of a previously con- 
jectured dynamical instability. 



Electronic bilayers exhibit a rich pattern of behavior, 
both on the static and on the dynamic level. While at 
high r s values (r s > r£ rvs ™) ) the bilayer is expected to 
crystallize (according to rf ystal > 20), and at very 
low r s values the RPA description is largely sufficient, the 
most interesting behavior occurs in the 1 < r s < r c a rystal 
domain, where the system is in the liquid state. This is 
the domain we focus on in this Letter. 

The most remarkable feature on the static level is 
that the system exhibits a series of abrupt structural 
changes as the ratio of the intcrlayer distance d and the 
2D Wigner-Seitz radius a is varied ||. These structural 
changes parallel the structural phase transitions in the 
solid phase but at finite temperature they are also 
combined with entropy increasing substitutional order- 
disorder transitions |3| where particles in layer 1(2) oc- 
cupy positions appropriate for particles in layer 2(1): this 
is signalled by the two pair correlation functions (PCF) 
hn(r) and h\2(r) becoming identical as d — > 0. 

The dynamical behavior, the collective mode structure 
in particular, has been studied by Swierkowski et al. [||, 
by Gold ||, by Zhang and Tzoar j7|, Golden, Kalman and 
collaborators ff-|l(| and by Moudgil et al. Jll| . (Some of 
these studies pertain to a superlattice [infinite number 
of layers], rather than to a bilayer; in qualitative terms 
the results can, however, be easily interpreted for the 
bilayer). Classical bilayer and multilayer structures that 
form in charged particle traps have also been studied the- 
oretically by Dubin [12a] and have recently been observed 
in ion traps [12b]. There are two problematic issues that 
make the predictions based on the calculations less than 
reliable |13| . The first issue relates to the approxima- 
tion technique used: most of the wor ks cited use 
methods which violate the 3 rd frequency-moment-sum- 
rule, whose satisfaction is well recognized [[t4| to be an 
important criterion for providing an acceptable descrip- 
tion of the collective mode behavior. The second issue 
concerns the use of the intralayer and interlayer correla- 
tion functions as inputs in all the calculations cited. No 
reliable pair correlation function (PCF) data - either for 
classical or quantum bilayer systems - have been avail- 
able until fairly recently: thus predictions of the collec- 
tive mode structure (which turns out to be extremely 
sensitive to the behavior of the inputted PCF) have been 
compromised from the outset. 

In this Letter we present the first consistent and re- 



liable calculation of the collective mode spectrum of a 
strongly coupled electronic bilayer liquid. Our results 
show features which are qualitatively different from the 
weakly coupled RPA results; they also show that earlier 
claims concerning t he p ossible emergence of a dynami- 
cal instability UlllllyJ cann °t be supported by a more 
consistent treatment of the correlations. The calculation 
is based on a purely classical model: (i) two 2D electron 
liquids separated by distance d; (ii) scattering on impuri- 
ties etc. neglected; (iii) no interlayer tunneling; (iv) the 
system is described as a binary liquid with interaction po- 
tentials ifu(r) = <p22 (r) = e 2 /r, <p\%(r) = e 2 / yj (r 2 + d 2 ). 
In addition, the model portrays a classical electron liq- 
uid where exchange and other quantum effects are ne- 
glected: this approximation is reasonable in the strong 
coupling domain where the particles are well localized. 
The intralayer coupling is characterized by the parame- 
ter T = e 2 /(aT) where a = l/y/(im) and T is the ki- 
netic energy per particle - the temperature in a classical 
system and (1/2)ef in a zero temperature 2D electron 
gas. Hence the equivalence T — > 2r s . The calculation 
of the dielectric matrix e^ I/ (kw) is carried out in the 
Quasi Localized Charge Approximation (QLCA), which 
has been applied successfully for the description of other 
strongly coupled Coulomb systems Jll]-[l7j] ; in the QLCA 
e(kw) becomes a functional of the intralayer and inter- 
layer PCF-s h\i{r) and /112(f) or of the corresponding 
structure functions S\\{k) and S\2{k): 

e(kuj) = I-uj 2 (ka)[uj 2 I-'D(k)}- 1 (1) 

(cJo = ( 27re " )i/ 2 ; the nominal plasma frequency of a 
single 2D layer) with |]| 

A7( k ) = 

^ £ gVbo-MStf (|k - q|) - % E MQ)Su(q)], (2) 
q I 

I is the identity matrix. In Cartesian space e(kw) and 
similarly D(k) are reducible to longitudinal e L (D ) and 
transverse e T (D T ) matrices. The dispersion relation for 
the longitudinal modes is then obtained from 

H^M^O (3) 

which leads to 

u? =cu ka(l±e~ kd )+ £>u(k) ±Df 2 (k) (4) 

With the neglect of retardation effects, the dispersion 
relation for the transverse modes is derived from 

|| £ T (k w )- 1 || = (5) 

which yields 

co 2 =D&(k)±£>f 2 (k) (6) 

The D-functions are to be expressed in terms of the struc- 
ture functions, according to Eq. |[ These latter have re- 
cently been calculated j| through the HNC (Hypernetted 
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Chain) integral equation for a wide range of T and d val- 
ues. Inputting them into the corresponding D-functions 
and using the latter in Eqs. [| and ^ one can generate 
a full description of the collective mode spectrum. The 
results are portrayed in Figs. 1-3 and the qualitative fea- 
tures of the collective mode dispersion are summarized 
below. (The figures are given for T = 40, corresponding 
to r s = 20: this r s value, while high enough for correla- 
tions to be dominant is within the domain of experimen- 
tal realizability.) 

1. The spectrum of collective modes comprises 4 
modes: 2 (longitudinal and transverse) in-phase modes 
(corresponding to the + sign in Eqs. | and |) and 2 
(longitudinal and transverse) out-of-phase modes (corre- 
sponding to the - sign in Eqs. | and |). 

2. The in-phase modes (where the two layers oscillate 
in unison) are not qualitatively different from the similar 
modes of an isolated 2D layer |ll|. In particular, for 
k — ► the longitudinal (plasmon) mode has the typical, 
quasi-acoustic uj ~ \f~k~ dispersion, while the transverse 
(shear) mode is acoustic, ui ~ k; both modes are softened 
by intralayer and interlayer correlations. 

3. The out-of-phase modes (where the oscillations of 
the two layers exhibit a 180° phase difference) are char- 
acterized by an u>(k = 0) > energy gap. The physical 
reason for the existence of an energy gap for layered sys- 
tems has already been discussed elsewhere [|| . The out- 
of-phase longitudinal mode in the RPA approximation 
has been identified as the acoustic plasmon j^,^] since 
for k — > 0, uj — > as k. The present calculation clearly 
shows the marked difference brought about by the strong 
correlations. Since at k = the isotropy of the system 
is unbroken, the plasmon and the shear modes share a 
common gap value. 

4. From Eqs. [| and ^| the gap value can be expressed 

as 



to 2 {0,d) 



d(qa)(qa) 2 e- qd S 12 {q). 



(7) 



With increasing d and consequently decreasing inter- 
layer correlations, w(0) shows a decreasing tendency and 
it virtually vanishes for d > 1.5 when the separated 
layers become practically uncorrelated J^]. This down- 
ward trend is, however, preceded by a slight upturn for 
< d < 0.16 (for T > 30). Although the details of 
this behavior are not well understood, it is most likely 
due to the substitutional disorder that prevails in this re- 
gion S: the eigenfrequencies of the localized modes in 
the substitutionally disordered phase are expected to be 
higher than in the substitutionally ordered phase p0[ . 
Within the domain investigated, the T-dependence of 
lu(0) is quite mild, but the QLCA being a strong coupling 
approximation, no inference concerning the behavior of 
u>(0) in the moderately coupled (r < 10) domain can be 
drawn from this observation. In fact, it is expected that 
for low enough T values, w(0) tends to zero to match at 
r = the predicted RPA behavior @. 



5. For finite k values all the four dispersion curves 
develop an oscillatory behavior, generated by the similar 
behavior of the inputted structure functions. This behav- 
ior has also been identified for the isolated 2D layer [16, 
17b]. The structure of the out-of-phase plasmon mode 
is of special interest here: the first sharp roton-like min- 
imum has attracted attention in earlier studies |^]|) 11 
which were based on the neglect or on a highly approxi- 
mate treatment of the interlayer correlations. It was sug- 
gested that the minimum of uj 2 may dip below oj 2 = [5a, 
9,11] or may, at least, reach the close vicinity of uj 2 = 
[5b, 5c]. The former behavior would indicate a dynami- 
cal instability (heralding the onset of CDW-type ground 
state), [5a, 11]; the latter has been interpreted as the on- 
set of a new high-fc, low frequency mode [5b, 5c]. Our 
results show that the roton minimum never drops below 
the value already reached by the dispersion curve of the 
2D layer. The consistent treatment of the interlayer and 
intralayer correlations thus precludes the existence of the 
effects conjectured in J5],[)],[ll] , virtually independently of 
any other approximation used. 

6. At high k- values, for a given d all the dispersion 
curves approach the same asymptotic frequency value, 
the frequency of a localized mode, a particle oscillating 
in the screening environment of the two layers: 

iv 2 (oo,d) = — j^2[fn(q)(Su(q) - 1) + V>i2{q)S 12 (q)]q 2 



d(qa)(qa) 2 [Su(q) - 1 + S 12 {q)e^ d ] (8) 



This result is, probably, only of academic interest, since 
it is unlikely that high-fc modes would survive the damp- 
ing mechanism operating in the system. 

The QLCA is not geared to describe damping processes 
and therefore our calculation fails to provide informa- 
tion on the damping of the collective modes. However, 
some qualitative statements can be made. We concen- 
trate on the out-of-phase modes only. There are three 
major damping mechanisms to be considered. These are 
(i) single pair excitation (Landau damping) , (ii) multiple 
pair excitations and (iii) diffusive-migrational damping 
p7[ . The effect of the single pair excitations can be eas- 
ily assessed from Figs. 1 and 2. Both of them show the 
pair excitation domain: it is clear that as long as the 
layer separation is not too large (d/a < 1.5) for small k 
values both the out-of-phase plasmon and shear modes 
are well outside the continuum and are thus immune to 
Landau-damping. In a highly correlated plasma multi- 
particle excitations are, however, also operative, with in- 
creasing importance at higher k values. The k — > oo 
portion of the dispersion curve emerging from the con- 
tinuum would probably be heavily damped by this pro- 
cess. The diffusive-migrational shifting of the quasi-sites 
p7| can originate from quasi-thermal diffusion or from 
tunneling between neighboring minima of the fluctuating 
potential. For T sufficiently high (r > 40, i.e. r s > 20, 
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probably close to crystallization the latter effect 

should be significantly diminished. The former, however, 
is sizable when its characteristic time becomes compa- 
rable with the period of the mode under consideration. 
Based on our earlier estimate [|l0| of the lowest surviv- 
ing oscillation frequency, u) m i n and on the MD results 
of Ref. (2^] (see also ]T^ ) concerning the lowest propa- 
gating wave-number value for the shear mode, one can 
conclude that for T = 40 (r s = 20) the d/a < 0.8 do- 
main and for T = 20 (r s = 10) the d/a < 0.6 domain can 
be safely assumend not to be seriously affected by this 
damping mechanism for either of the modes, (cf Fig. 
1/a). For higher d values the gap frequency u>(0) is be- 
low (jJmin and thus the shear mode would be propagating 
for lo > LUmin only; in the domain u) < Lo m in the longitu- 
dinal mode would be stripped of its correlational features 
and would revert to an RPA type acoustic plasmon with 
u(k -> 0) -> 0. (cf. Fig. 1/b). 

The collective mode structure of the strongly coupled 
bilayer liquid presented in this Letter bears a close re- 
lationship to the phonon spectrum of the bilayer solid, 
recently calculated by Goldoni and Peeters 0]. The 
four modes in the solid phase can be identified as the 
transverse acoustic and the longitudinal quasi-acoustic 
(~ y/k) phonons and the transverse and longitudinal op- 
tical phonons. The "gaps" exhibited by the latter are, 
in general, different because of the anisotropy of the lat- 
tice. In contrast, in the liquid state, there is only one 
single isotropic gap, as determined in this Letter. This 
liquid gap value is typically slightly above the arithmetic 
average of the optical frequencies of the transverse and 
longitudinal phonons, in the solid phase. 

Concerning the possible observation of the features 
predicted in this Letter, one can suggest three main ar- 
eas that should lend themselves to direct experimental 
verification: (1) the existence and the nonmonotonic in- 
dependence of the k = energy gap; (2) the existence 
of a transverse shear excitation with a high frequency 
and expected low damping (this is in a sharp contrast 
to the usual scenario for the shear mode in the liquid 
phase, which vanishes for k — » |l^ , |2~i| , |22]| ; (3) the non- 
existence of the predicted PJ^j9|,pl|instability or low fre- 
quency mode in the vicinity of the first roton minimum. 
We note that the reported Raman scattering experiments 
are inconclusive because of the low r s and relatively 
high k values involved. Recent advances in fabricating 
high r s samples [|4j and small layer separation should 
render the suggested experiments feasible. 

In summary, we have obtained a comprehensive picture 
of the collective mode structure of an electronic bilayer in 
the strongly coupled liquid phase. This structure is qual- 
itatively different both from that of the weakly coupled 
bilayer electron gas (describable in the RPA) and from 
the phonon spectrum of the bilayer solid and exhibit a 
number of remarkable features which should be exper- 
imentally verifiable. The calculation avoids the pitfalls 
of earlier approaches which stem from the inconsisten- 
cies in the approximations used for the interlayer and for 



the intralayer PCF-s, and from the violation of the 3 
frequency-moment-sum rule. 
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FIG. 1. The four principal modes for T = 40 (r s = 20): 
(a) d/a — 0.3 (b) d/a — 1.0. The shaded region is the pair 
excitation continuum. For d/a = 1.0 (ui m in ~ 0.45o;o) the dif- 
fusional domain in which the shear mode is non-propagating 
and the out-of-phase plasmon exhibits an RPA-like behavior 
is indicated. 

FIG. 2. The gapped out-of-phase ("acoustic") plasmon 
for different layer separations. The shaded region is the 
pair excitation continuum. In the u < u> min domain the 
d/a = 1.0, 1.5, 2.0 gap values are spurious and the dispersion 
should resemble the RPA acoustic type behavior. 

FIG. 3. Gap (w(0)) value for r = 40 and T = 30 as func- 
tions of the layer separation. Below the estimeted 0J m i„(r) 
value the gap is spurious. The inset compares w(0) and w(oo) 
for T = 40. 
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